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Abstract 

o ■ 

i— I I Recently, the contributions of chiral logarithms predicted by quenched chiral per- 

turbation theory have been extracted from lattice calculations of hadron masses. 
. We argue that a detailed comparison of random matrix theory and lattice calcula- 

tions allows for a precise determination of such corrections. We estimate the relative 
size of the mlog(m), m, and m 2 corrections to the chiral condensate for quenched 
^ SU(2). 
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The identification of logarithmic corrections in the quark mass predicted by 
quenched chiral perturbation theory [1,2] in lattice gauge results is a long 
standing problem. It seems that the latest numerical results [3-6] on hadron 
masses in quenched lattice simulations allow for an approximate determination 
of these log(m) contributions. The determination of these logarithms is an 
important test of chiral perturbation theory which in turn plays a central role 
for the connection of low-energy hadron theory on one side and perturbative 
and lattice QCD on the other. 

In a completely independent development, it has been shown by several au- 
thors that chiral random matrix theory (chRMT) is able to reproduce quan- 
titatively microscopic spectral properties of the Dirac operator obtained from 
QCD lattice data (see the reviews [7,8] and Refs. [9-12]). Moreover, the limit 
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up to which the microscopic spectral correlations can be described by random 
matrix theory (the analogue of the "Thouless energy") was analyzed theo- 
retically in [13,14] and identified for quenched SU(2) lattice calculations in 
[15]. 

The following analysis uses the scalar susceptibilities, so we first give their 
definitions. The disconnected susceptibility is defined on the lattice by 

<*. L/f 1 \ ' \' (1) 

Xlattic N \ 2^ (iAt + m)(Uj + m) j N y A, + m J , V J 



where N = L A denotes the number of lattice points and the are the Dirac 
eigenvalues. After rescaling the susceptibility by NT, 2 (E = absolute value of 
the chiral condensate for infinite volume and vanishing mass) chRMT predicts 
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where the rescaled mass parameter u is given by u = mTL A . (For details we 
refer to [15].) 

We shall also use the connected susceptibility which is defined on the lattice 
by 
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Fig. 1. The ratio of Eq. (6) for the scaled susceptibilities plotted versus mSL 2 (in 
lattice units) for (3 = 2.0 and four different lattice sizes, N = 4 4 , 6 4 , 8 4 , and 10 4 . 

The chRMT result reads 

^ = AuK^u) J ds (1 - s)I (2su) . (5) 



Fig. 1 presents the deviation of the (parameter-free) random matrix prediction 
from the lattice result, more precisely the ratio 

ratio = (xi a ttice - Xrmt) / (xrmt) , (6) 



where x can either be the disconnected (only this choice was investigated in 
[15]) or the connected susceptibility. 

The motivation for investigating ratio rather than xiattice itself is that in 
Eq. (6) finite size corrections cancel to a remarkable degree, allowing us to use 
data from smaller m values. We have seen in Fig. 2 of [16] that the knowledge of 
finite size effects which we gain from RMT allows us to find the thermodynamic 
limit of the chiral condensate from extremely small lattices. This can also be 
formulated in the following way: for a given value of ratio in Fig. 1, the finite 
size corrections for all four lattice sizes are expected to be similar, as the 
corresponding values of m^L 2 oc mL 2 are very close, which is why we have 
plotted ratio against mSI 2 in Fig. I. 

What do we expect beyond the Thouless energy? Then, the lattice is large 
enough so that the valence pion, which is the lightest particle, fits on the 
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lattice. Naturally, all other particles also fit on the lattice, and therefore we 
expect that the chiral condensate and the two susceptibilities will rapidly 
approach their thermodynamic limit. 

For a finite lattice and a non-vanishing mass, the chiral condensate is given 
by 



N 



^lattice (m) 



N \~ i\ k + m i 



(7) 



In the quenched theory, the connected susceptibility is given simply by 



(8) 



so we can find the infinite-volume behavior of ^ conn from that of a. We expect 
from chiral perturbation theory [17] that the chiral condensate has the form 



a(m) = £ 



1 — Am log(m) + Bm + -Cm 2 + 



(9) 



Eq. (9) requires several comments. In the continuum, quenched chiral per- 
turbation theory predicts a leading term proportional to ^r-log(m), where 
(z/ 2 )/L 4 is the topological susceptibility [17, Sec. 7]. We argue that this lead- 
ing term should be absent in our case. For finite lattice spacing the Atiyah- 
Singer-index theorem does not apply for staggered fermions. Therefore the 
role of topology has to be interpreted with care. We have seen in [11] that 
the small Dirac eigenvalues are well described by random matrix results for 
v = 0. This means that the quasi- zero modes related to topology are shifted 
to such large values that they are not visible. (This is presumably due to 
discretization errors proportional to a 2 , with a the physical lattice spacing.) 
Thus the violation of axial symmetry which generates the logarithmic term in 
the quenched case is dominated by the explicit quark masses, which motivates 
Eq. (9). It would be very interesting to study the v ^ sector for which we 
expect a leading log(m) term, which might require, however, very small a and 
a large number of lattice points. 

Eq. (9) implies that in the thermodynamic limit 

XiattL = S [-A log(m) - A + B + Cm + • • •] . (10) 



On the other hand, the large-volume limit of the RMT susceptibility is 

conn ' _ 

XRMT - 4m2L 4 • 
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Putting the two expressions together, we find that 

4 

ratio -> (mEL 2 ) 2 - [-A log(m) - A + 5 + Cm+ • • •] - 1 . (12) 

2j 



Strictly speaking, the —1 ought to be neglected in comparison with the first 
term as L — > oo. However, we should be prepared to observe some sub- leading 
corrections in the data taken on finite lattices. 

We have confronted Eq. (12) with lattice Monte Carlo data for two values of 
the coupling strength /3, f3 = 2.0 and f3 = 2.2. The lattice sizes and numbers 
of configurations are given in Table 1. 

Table 1 

Lattice sizes and numbers of configurations for the lattice data. 





P = 


2.0 






L 

# of configs 


4 

49978 


6 

24976 


8 

14290 


10 

4404 




= 


2.2 






L 

# of configs 


6 

22292 


8 

13975 


10 
2950 


12 
1388 



To check Eq. (12) we did the following for both values of f3: 
We chose different values for ratio = foj and determined the values of mEL 2 
for which they were reached for our different lattice sizes. Let us denote these 
numbers by Y(L, Eq. (12) implies that 



Y(L,bi)< 
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where r(b) will be proportional to 1/b as b 
large values of b, we used the ansatz 
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to fit our data. In Eq. (14) not only Y 2 has statistical errors, but also m. In 
our x 2 fit) however, only the errors of Y~ 2 are taken into account. 

Obviously, the values Y(L, bi) for the same lattice size L are highly correlated. 
It is, however, unclear how to calculate the correlations of these quantities, 
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Fig. 2. The value of Y = mUL 2 at (5 = 2.0 for which ratio = b for various values 
of b as a function of m. Larger values of b belong to smaller values of 1/Y 2 . The 
rightmost filled dots correspond to L = 6, the leftmost to L = 10, whereas the open 
dots represent data for L = 4, which were not used in the fit. All quantities are 
measured in lattice units. 

which are related to the original lattice results only in a rather implicit manner. 
Moreover, correlated fits tend to have problems [18]. Therefore we decided to 
ignore correlations completely, although this will lead to an underestimation 
of the errors on the fit parameters. 

For the thermodynamic limit of the disconnected susceptibility we assume the 
same form as Eq. (10). In RMT, the large- volume limit is given by 

X ™ T 8m^ (15) 



so that the ansatz of Eq. (14) applies as well. 

In Figs. 2 and 3 we plot Y~ 2 versus m together with the fits for (3 = 2.0 
and 2.2, respectively. In the case of the connected susceptibility we used 6, = 
2.0, 3.0, 4.0, 5.0 (P = 2.0) and k = 5.0, 6.0, 7.0, 8.0 {(5 = 2.2) and obtained the 
results of Table 2. For the disconnected susceptibility we used bi = 1.0, 2.0, 3.0 
(P = 2.0) and h = 6.0, 7.0, 8.0 (f3 = 2.2) and found the values given in Table 3. 



The main message of Figs. 2 and 3 is that without any doubt the data are 
not fitted by horizontal lines. This demonstrates the presence of additional 
contributions in the quark mass. The approximate linearity of the curves for 
small m shows that the logarithmic contribution is the dominant one. For the 
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Fig. 3. Same as Fig. 2 but for (3 = 2.2 and with the dots, from left to right, 
corresponding to L = 12, 10, 8, and 6, respectively. 

connected susceptibility, the data are well fitted by the ansatz (9), i.e., with 
only the three leading corrections. For the disconnected susceptibility, our 
statistical precision does not allow for a precise determination of the ratios 
B/A and C/A. For very small lattices (4 4 in Fig. 2) finite size effects seem to 
spoil our analysis. 

It is clear from Figs. 2 and 3 that one would really like to have numerical 
simulations with substantially larger statistics and larger lattices. As the ap- 
plicability of RMT to the description of the low-energy Dirac spectrum is by 
now well established we can limit ourselves in the future to the calculation 

Table 2 

Fit parameters for the connected susceptibility. 



p B/A C/A q s X 2 /dof 

2.0 2.29 ±0.63 -5.97 ± 5.17 43.9 ± 4.5 0.25 ± 0.02 0.50 

2.2 0.86 ±0.18 -2.46 ± 1.86 486 ±19 0.81 ± 0.05 1.75 

Table 3 

Fit parameters for the disconnected susceptibility. 



(3 B/A C/A q s X 2 /dof 

2.0 1.9 ±3.1 -12 ±32 31 ± 16 0.05 ± 0.05 0.02 
2.2 -1.45 ±0.48 18.7 ± 4.1 569 ± 127 -0.60 ± 0.61 0.28 
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of just the lowest eigenvalues instead of the complete spectrum. This should 
allow us to gain the necessary statistics. 

To conclude, let us remark that the aim of this paper is primarily to draw 
attention to this new method to extract chiral logarithms and other corrections 
in the quark mass, and to stimulate the discussion of their interpretation. The 
obvious next step is to analyze the susceptibilities within the framework of 
quenched chiral perturbation theory. 
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